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Abstract 

A perturbation method is presented which can be apphed to the description of a wide range 
of physical problems that deal with dynamics of dipolarly coupled spins in solids. The method 
is based on expansion e-^"'"'^ in a perturbation series. As example of application the method, the 
multile-quantum coherence dynamics in three and four spin cluster are considered. The calculated 
OQ- and 2Q intensities vs the duration of the preparation period give closed agreement with exact 
results and simulations data. The exact solutions for Jqq and J2Q coherences in four spin system 
are obtained. 
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Dipolar coupling spin dynamics represents a significant interest from both the point of 
view of the general problems of physics of many bodies, and from a position of the nuclear 
magnetic resonance (NMR) ^,3]. In solids, the evolution of a spin system under the 
dipole-dipole interaction (DDI) involves many spins and leads to unsolvable problems. Even 
the analysis using the numerical calculation becomes difficult because the number of states 

= 2" is growing exponentially with the increasing of n. So in existing theories, only 
macroscopic characteristics such as spin-spin relaxation times, the second and the fourth 
moments of resonance lines were taken into account [l^. This difficulties are very clearly 
displayed in multiple-quantum (MQ) spin dynamics. The MQ phenomena involves various 
multiple-spin transitions between the Zeeman energy levels and form MQ coherence at times 
t > uj^^ , where ujd is the characteristic frequency of DDI P]. The problem in analytical 
description of MQ is that the different modes of coherence grow in at different times, with 
higher modes requiring longer excitation times than lower modes ^. Hence, uJdt > 1 is 
not a small parameter, and, at first glance, perturbation theory methods cannot be used 
to study MQ dynamics. Indeed only simple exactly solvable models of spin system such 
as two and three dipolar coupling spin-1/2 0,0] or one- dimensional linear chains spin [3] 
system were analyzed theoretically. The last achievement in this direction is the model with 
identical DDI coupling constant for all spin pairsj^, |^. Note, that the simplified calculations 
are essential for the case of indentical DDI coupling constant has been already mentioned 
[lfll |. Such approaches cannot describe MQ processes except for zeroth and second-order 
coherences. Thus, the development of method that can successfully represent important 
features of MQ dynamics with a larger then OQ and 2Q coherences is needed. 

Importance of the analytical description of the MQ processes in solids is that involving the 
excitation of collections of dipolar coupling spins, can provide important structural, as well as 
the spin dynamics informations. Moreover, during the past few years, the NMR is considered 
as the best candidate 12 , ^ for experimental realization quantum information processes. 
It was experimentally demonstrated the creating pseudopure spin states in large clusters 
of coupled spins by MQ method and that dynamics of the quantum entanglement is 
uniquely determined by the time evolution of MQ coherences 6]. Thus, investigation of the 
quantum information processes with MQ methods are of current interest. 

We present a perturbation approach to the problem of dipolar coupling spin dynamics 
in solids. DDI of all spins was divided into several groups that are characterized by the 
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identical DDI coupling constants. Since the magnitudes of the dipolar coupling constants 
vary inversely with the cube of internuclear distance, the coupling constants are different 
for these groups. 

Our main idea is to take into account in MQ NMR dynamics influence of the groups 
with different degree of accuracy. Spin groups with smaller DDI coupling constants can 
be considered as perturbation (spins located far apart) while the nearest neighbours are 
taken into account exactly. As the result, we can develop a perturbation method that allows 
obtaining the description of the MQ with a large coherence evolution under DDI in an 
analytical form. The proposed approach will be a powerful method to describe wide ranges 
of physical problems that deal with dynamics of dipolar coupled spins in solid. On the one 
hand, this approach uses the advantages of exactly solvable models 0, . On the other 
hand, it simplifies calculations by using a perturbation technique. Results will show that the 
perturbation method can be applied to solve complex spin-dynamics problem and to obtain 
the solution in an analytical form. The method is based on the differential method Q| 
expresses e'^^'^ as an infinite product of exponential operators ^|. In the case when the 
norm of operator B is small then one of operator A, \\B\\ < \\A\\, we will try to obtain the 
perturbation series which takes into account up to second order terms in ratio / ||^||. 
Then the problems in description the MQ dynamics will be significantly simplified. 

Let us consider a spin system with Hamiltonian, 7i which includes only two parts with 
different DDI constants a and (3: H = A + B , where A = aA and B = j3B, a = \\A\\ and 
P = \\B\\ are the norms of operators A and B, respectively {a > f3 and [^, -B] 7^ 0). The 
evolution of the spin system is governed by propagator 

We seek to express (1) as a series in such that 

^-UiaA+m =e-itl3B^^^^^ (2) 

where operator a (t) obeys the differential equation 
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= (^) ^ (^) ^ (3) 

with initial condition 

a (0) = 1 (4) 
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and (t) = e-'*f^^ Ae'*^^ . Assume that at > 1 and f3t < 1. First, we will restrict ourself 
by keeping only first-order terms that are linearly proportion to /3. This leads to 

i^^^^^^ = a (A - itf3 [B, A]) (7(°) (t) . (5) 

To solve Eq.(5) we will use the iterative method. We will search the solution of Eq.(5) as a 
series on parameter (at)"' : 

oo 

aA{t) = J2a'^^\t) (6) 

n=0 

Taking into account that cr^^ (t) — 1 the solution of Eq.(5) for n = 1 can be obtained: 

(t) = (-iat) A. (7) 



For n — 2, ones obtains 



"Tm-^-^f^lA' + ^lB.A]] (8) 



Keeping only linear in ^ terms the following expression for operator o"^^ (t) can be obtained 



^ ( A'" - ^ I (m - 1) - J] A^'-^-^BA^ ) ) (9) 



m-2 



Using Eq. (9) we obtain 



^-it(aA+m ^ ( 1 _ ^ y a"7^^^ y 7/^ ^"{[ a^ba^-^ 1 e 



(10) 

where [B,A]j^_^^ denotes the repeated commutators [[[... [[S, 74],74] ...A]^_^-^ . After summa- 
tion over j in Eq. (10) we have 



\ n=0 ^ ' J 

where 

{S,^°} = 5and = (12) 

In the limit as the number of steps A?" — > oo we obtain that limjv=oo ^(jv+2)! ) ~ ■ 
Consequently, the exponent in (11) can be presented in the limit as A^ — > oo in the following 
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-i( QfA+g '""V'^r [B,A]n+i I _■ tA 

form: lim7v=oo e ^ " ^ ^ / = e , which does not include any terms with 

B, and the summing over n up to indefinite, results in 

^-^tiaA+f3B) = (^1 - j^^ dxe-"^5e"^^ C"*"*^, (13) 

which is a well known formula for expansion an exponential operator in a perturbation 

series To obtain expansion containing only linear to - terms we have to require that 

1 

a ^lN+2y. '^^^^ requirement imposes restrictions also on time: t ^ (^^ {N + 2)!^ 

1 

So for the smallest of times t <^ a + 2)!^ or t <^ the Eq. (10) includes only 

terms linear m j3 . In an analog way, we obtain the expansion up to second order in the 



ratio - : 

a 

^ a '^^m\k\(k + m + l) 

m=0 fc=0 ^ ^ 

(-l)'"+^(zx)^+^+"^+P+^ A'+'^EAPEyl^ 
^lai m!A;!/!p!(/+p+l) (A; + / + m +p + 2)^-^^ ^ ^ 

The series expansion (14) can be used not only for the small parameter ^ < 1, but 
independently, for the parameter x = at. Formula (14) can be easily generalized for a case 
when the exponential operator contains arbitrary number of the non-commutative operators 
and can be extended to include various power of the operators. Eqs. (14) appears to be 
complex at glance, but in fact it is quite simple to use, as the following examples will 
illustrate. 

Let us consider a cluster of three dipolar-coupled spin-| nuclei. The MQ dynamics in 
the rotating frame is described by propagator (1), where the time-independent average 
Hamiltonian is given by 

j<k 

and Ij' and I J' are the raising and lowering operators for spin j. The dipolar coupling 
constant, djk, for any pair of nuclei j and k in the cluster, is given by 

= (1-3 COS 9jk) , (16) 

jk 

where 7 is the gyromagnetic ratio of the nuclei, rj^ is the internuclear spacing, and 6jk is 
the angle the vector fjk makes with the external magnetic field. In the high-temperature 



approximation the density matrix at the end of the preparation period is given by 

= e-'«V(0)e^^* (17) 
where p (0) is the initial density matrix in the high-temperature approximation 

P(0) = E^I' (18) 
i=i 

/J is the projection of the angular momentum operator on the direction of the external field 
for an spin j. The average Hamiltonian (15) can be divided into the three parts according 
to the number of the different coupling constant di2 > d23 > c^is '■ 

H^Hu + H2s + H,-i, (19) 

where 

Hjk = + lilk) with J ^ k and J, k^l,2, 3. (20) 

The experimentally observed values are the intensities, J„q (t) of multiple-quantum coher- 
ences: 

-^riQ it) = y^^2 (0) Yl pIi n = m^p- m,g, (21) 

where m^p and rrizq are the eigenvalue of the initial density matrix (18). The perturbation 
method described above is used to calculate the time evolution of MQ coherences. Using 
expansion (14) with aA — H12 and j3B — H23 + i?i3 and keeping terms up to eighth order 
in X — at , the normalized 0-quantum ( Jqq) and 2-quantum (J2q) intensities are given by 

32x^ 256x^ 512x^ 
^OQ - 1 - ^ + —p, T—- + 



3 9 125 945 
/5\V8a;2 40a;^ 176x^ 15Ux^ 



and 



Ax^ l&x^ 128x6 256x^ 

= — + — 7^ + 



3 9 135 945 

(Ax^ 32a;^ 128a;*^ 1024a; 



where (^)^ = + • Let us compare formulas (22) and (23) with results from Eq. 

(14) which the terms with x — at will be taken into account exactly. By summing over n, 



m, /, and k up to indefinitely in (14), we obtain the analytical expressions of the intensities 
of - quantum 



JoQ = cos 4x - 2 



and of 2- quantum 

sin 2x 



J- 



2Q 



2x 



+ 2 ( — ) (cos X + cos 3x — cos Ax — x sin Ax) ^ (24) 



/3 



a 



cos 2x + 2 



[3 



a 



/3 



a 



cosx 



smx 



(25) 



Now let us compare intensities (22) - (25) with the exact solution [5|. Figs. 1 and 2 show 
the evolution of the normalized OQ and 2Q coherences for three spin cluster, where ^ = 0.3 
and at t = the spin system is in thermal equilibrium (18). All approaches, perturbations 
(Eqs. (22) -(25)) and exact give closed agreement up to x = 0.75 (in unit of ^), both for 
OQ - and 2Q- coherences. The exact account of influence of the nearest neighbours gives a 
good agreements up to x = 2. 

As a second example we consider is cluster consists of four spin arrange in corners of 
square in an external magnetic field penpendicular to the square plane. In this case the MQ 
spin dynamics is described by the average Hamiltonian 



H = Hi + H2, 



(26) 



with two different dipolar coupling constants Di and D2 , where Di and D2 are the dipolar 
coupling constants between nearest neighbors and spins at opposite sites, respectively (£ = 



and 



^1 = ("t) ^ + ^^^^^ = 

V / j=i 

V / j=i 



(27) 



(28) 



Using the expansion (14) up to eighth order in x = at, the normalized 0-quantum (Jqq) 
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Joo = 1 - + -x'^ - —x^ + —x^ 
OQ 4 18 28 

and 2- quantum (J2Q) 



PY f^^2 13 4 121 6 



a 



X X + 

6 60 



X 



—x' ] (29) 
630 ' ^ ^ 
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(30) 



7 



intensities can be determined. Formulas (29) and (30) will be compare with results from 
Eq.(14) in which the terms describing interaction of the neighbour spins will be taken into 
account exactly. By summation over n, m, I, and k up to indefinitely in Eq.(14) we obtain 
the analytical expressions of the intensities of - quantum 



^ (^l + sin2 2a; + 2sin2y2a;) + y 1^^^ ^2 sin^ 2a; - ^ sin 2722; j (31) 



and of 2- quantum 
J2Q 



~ (^sin^ 2x + 2 sin^ V2x^ + ^ (^^^ ^2 sin^ - ^ sin 2\/2a;j (32) 
coherences. 

To control the perturbation results (29) - (32) we obtained the exact solution for Jqq 



JoQ = ^ - icos |^2a;y'2 + ^-^ j + ^cos4xcos (^2a;-^ - cos^ ^x-^ (33) 
and for J2Q 




= -7 + 7C0S I 2x'\/2+ { ^ ) I + -sin^2xcos (^2a;^J + -cos^ (^x^J (34) 

coherences and fulfilled the numerical analysis of the MQ dynamics. The exact solutions 
(33) and (34) and computer simulation of the MQ coherences of four spins cluster have been 
obtained with a PC using the MATLAB package. 

Figs. 3 and 4 show, that perturbation results (29) and (30) are in the agrement with 
(31) and (32) and the exact solutions (33) and (34) and with the simulation data one up to 
X = 1 (in unit of ^). Calculations in which the interaction between the nearest neighbours 
is taken into account exactly (Eqs (31) and (32)) are in close agreement with the exact 
solutions (Eqs. (33) and (34)) and simulation data up to a; = 3. 

In conclusion, a perturbation method was developed which is based on the expansion of 
operator exponent in a perturbation series. Then the perturbation approach was applied to 
the description the MQ spin dynamics in sohds. The analytical expressions for OQ and 2Q 
dynamics in a three and four spin clusters in solids were obtained. In the four spin cluster 
the exact solution was obtained. The calculated OQ- and 2Q intensities vs the duration of 



the preparation period agree well with exact solutions for three j5| and for four spin clusters 
(Eqs.(33) and (34)). 

The developed method can be extended to include various power of the operators with 
small norm and applied to the description widely range of physical problems deal with 
dynamics of dipolar coupling spins in solids. The results in an analytical form can be 
use to extract from experimental data the dipolar constants and the molecular structure 
information. 
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Captions for figures. 
Figure 1 

Fig.l Time dependences (in units of ^ ) of the normalized intensities of OQ coherence 
(solid-fine is exact solution [? ], dot-fine is the calculation using Eq.(22) and dash-fine is the 

calculation using Eq.(24)) 
Figure 2 

Fig.2 Time dependences (in units of ^ ) of the normalized intensities of 2Q coherence 
(solid-fine is exact solution[? ], dot-line is the calculation using Eq.(23) and dash-fine is the 
calculation using Eq.(25)) 

Figure 3 

Fig. 3 Time dependences (in units of ^ ) of the normalized intensities of OQ coherence 
in four spin cluster. Solid-fine is exact solution (Eq.(33), dot-line is the calculation using 
Eq.(29), dash-line is the calculation using Eq.(31), and open circle is computer simulations. 

Figure 4 

Fig. 4 Time dependences (in units of ^ ) of the normalized intensities of 2Q coherence 
in four spin cluster. Solid- line is exact solution (Eq.(34), dot-line is the calculation using 
Eq.(30), dash-line is the calculation using Eq.(32), and open circle is computer simulations. 



10 



1.0-, 



0.8- 




0.0 0.5 1.0 1.5 2.0 



X=at 




X=at 




X=at 



